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The magnetic field opens a gap in the edge state spectrum of two-dimensional topological insula¬ 
tors thereby destroying protection of these states against backscattering. To relate properties of this 
gap to parameters of the system and to study dynamics of electrons in edge states in the presence 
of inhomogeneous potentials, the effective Hamiltonian theory is developed. Using this analytical 
theory, quantum-mechanical problems of edge-state electron transmission through potential steps 
and barriers, and of motion in constant electric field are considered. The influence of magnetic 
field on the resistance of two-dimensional topological insulators based on HgTe quantum wells is 
discussed together with comparison to experimental data. 

PACS numbers: 73.21.-b, 73.23.Ad 


I. INTRODUCTION 

The two-dimensional (2D) topological insulators (also 
known as quantum spin Hall states) experimentally re¬ 
alized in HgTe-based quantum wells (THU represent a 
unique state of matter characterized by the existence 
of a pair of helical counter-propagating one-dimensional 
(ID) conducting channels at the edges of 2D system. 
These edge states have gapless energy spectrum and are 
topologically protected due to time reversal symmetry, 
which means that electrons in these states move with¬ 
out backscattering if the inelastic scattering is neglected. 
When the Fermi energy stays in the bulk insulating gap, 
the edge states determine electron transport. The resis¬ 
tance quantization observed 0 at low temperatures in 
the HgTe quantum well samples of macroscopic (1 /iin) 
length and subsequently confirmed in nonlocal resistance 
measurements 0 is a convincing proof for ballistic trans¬ 
port of electrons in the edge states. 

The problem of influence of magnetic fields on the 
transport properties of 2D topological insulators has ap¬ 
peared immediately after the discovery of these systems. 
Experimentally, it was found [l|-0 that the field perpen¬ 
dicular to the quantum well plane, B z , considerably sup¬ 
presses the conductance already at B z ~ 0.01 T, while 
the in-plane magnetic field suppresses the conductance 
only when it is strong enough, of the order of 1 T. It was 
also established that the conductance depends on the di¬ 
rection of the magnetic field in the quantum well plane 
0. Later on, it was demonstrated @, 0 that similar 
phenomena also take place in long samples, where the 
edge-state transport is not ballistic and the resistance is 
much larger than the resistance quantum. However, the 
effect of perpendicular magnetic field in these conditions 
is much smaller and manifests itself at B z ~ 1 T, sim¬ 
ilar to the case of in-plane field. To describe the effect 
of the perpendicular field, several theoretical approaches 
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have been proposed, including numerical simulation of 
transport based on a 2D disordered site model 0, as 
well as analytical consideration of the possibility of en¬ 
hanced backscattering due to edge spectrum nonlinear¬ 
ity induced by the magnetic field 0 and of interference 
effects at a disordered edge allowing for loops of the he¬ 
lical edge states 0- The approach of Ref. 8 gives an 
adequate description of the conductance suppression phe¬ 
nomenon under the assumption of strong disorder, when 
the amplitude of random potential perturbation exceeds 
the bulk gap energy. The behavior of the conductance 
under in-plane field has not been studied in detail, al¬ 
though numerical simulations 0 also show a suppression 
of the conductance under the action of this field. 

A convenient analytical method for theoretical stud¬ 
ies of edge-state transport in 2D topological insulators 
is based on the effective Hamiltonian theory, which as¬ 
sumes reduction of the 2D quantum-mechanical problem, 
which essentially requires numerical solution in the pres¬ 
ence of disorder potential, to a ID problem of electron 
motion along the edge. The effective 2x2 Hamiltonian 
for edge states 11, [l2[ has the form Hvka z , where v is 
the edge-state velocity, k is the operator of electron mo¬ 
mentum along the edge and a z is the Pauli matrix. This 
Hamiltonian describes two branches of edge states with 
gapless energy spectrum. However, it does not take into 
account effects of magnetic field and potential perturba¬ 
tions. Meanwhile, it is well understood 0 m that exter¬ 
nal magnetic field (or the presence of magnetic impurities 
in the system) violates the time reversal symmetry and 
makes the edge states no longer gapless. The gap in the 
edge state spectrum means that the states with opposite 
spin projections mix with each other and, therefore, elec¬ 
tron scattering between counter-propagating states be¬ 
comes possible. Such a backscattering should reduce the 
conductance in a similar way as it takes place in con¬ 
strictions [HI or finite-width strips of 2D topological in¬ 
sulators, where the gap is formed because of the mixing 
of states from the opposite edges due to overlap of their 
wave functions [lij . Therefore, inclusion of the gap open¬ 
ing effect of magnetic field into the Hamiltonian of edge 
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states is a necessary step in developing the analytical 
theory of magnetotransport in 2D topological insulators. 
The most general form of such a 2 x 2 matrix Hamil¬ 
tonian has been proposed in Ref. 12. It differs from 
the zero-field case by the presence of an additional term 
y gij aijBi&j which is linear in components of magnetic 
field Bi (here i and j are Cartesian coordinate indices). 
For application of this Hamiltonian to HgTe quantum 
wells, the coefficients a t j , which determine the effect of 
the magnetic field, have been calculated numerically Hi¬ 
lt is desirable, however, to have a more convenient rep¬ 
resentation of the effective Hamiltonian, where these co¬ 
efficients are expressed through the known parameters of 
the HgTe quantum well system. 


The basic goal of the research given in this paper is 
to develop the all-analytical effective Hamiltonian the¬ 
ory for ID motion of electrons in the edge states in the 
presence of magnetic field, spin-orbit interaction, and po¬ 
tential perturbations. The key point of the forthcoming 
consideration is a transition from the two-subband 2D 
Hamiltonian of electrons in HgTe quantum wells EE to 
an effective ID Hamiltonian for a couple of edge states. 
The influence of magnetic field is taken into account by 
using perturbation theory. The equation (7) derived in 
this way is similar to the massive Dirac equation de¬ 
scribing ID motion of spinless particles with a two-band 
energy spectrum. The gap between the bands in this 
spectrum considerably depends on the direction of the 
magnetic field and vanishes for a special orientation of 
the field. To show how the presence of the gap influ¬ 
ences ballistic propagation of electrons in the edge states 
and how the backscattering probability depends on this 
field, a straightforward application of the theory to stan¬ 
dard quantum-mechanical problems, such as transmis¬ 
sion through potential barriers and motion in a uniform 
electric field, is carried out. Edge-state transport in 
the random potentials is also discussed. It is demon¬ 
strated that the gap in the edge state spectrum leads 
to a quadratic dependence of the resistance on the mag¬ 
netic field strength. Some experimental magnetotrans¬ 
port data for 2D topological insulators, in particular, ini¬ 
tial quadratic increase of the resistance as a function of 
the in-plane magnetic field and dependence of the resis¬ 
tance on the orientation of this field are explained within 
the proposed theory. 


The paper is organized as follows. Section 2 is de¬ 
voted to derivation of the equation for ID motion in the 
edge states. Application of this equation to calculation 
of transmission and reflection probabilities, necessary for 
evaluation of the ballistic conductance via edge states, 
is done in Sec. 3. A discussion relating the results to 
experimental data on the conductance of 2D topological 
insulators based on HgTe quantum wells is presented in 
Sec. 4. 


II. HAMILTONIAN FOR EDGE STATES 


The effective 4x4 matrix Hamiltonian describing a 
2D insulator phase in symmetric HgTe quantum well in 
the basis of two subbands (one of them is the interface¬ 
like subband formed as a result of hybridization of 
conduction-band states with light hole states and the 
other is the ground heavy hole subband, with quanti¬ 
zation energies E\ and £2 respectively) was presented by 
Bernevig, Hughes, and Zhang [hj and improved in subse¬ 
quent works EI EM to account for spin-orbit cou¬ 
pling and magnetic fields. This Hamiltonian is written 
below in the approximation neglecting inessential small 
diagonal terms quadratic in electron momentum: 
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the quantum well plane, ip r is the electrostatic potential, 
A r = £1 — £2 is the gap energy whose possible variation in 
the plane is described as A r = A + <5 r , where A is the av¬ 
eraged gap, k± = k x ±ik y = — id x — eB z y/Kc±dy , e is the 
absolute value of electron charge, and B± = B x ±iB y . To 
describe the orbital effect of the perpendicular magnetic 
field B z , the vector potential is chosen as (— B z y , 0). The 
coefficients A, K, K z . and D are parameters of the ef¬ 
fective 2D Hamiltonian, they can be obtained directly by 
applying the Kane Hamiltonian for calculation of eigen¬ 
states in the quantum well system (see Ref. 16 and ref¬ 
erences therein). In particular, A describes coupling be¬ 
tween the subbands at finite 2D momentum, K and K z 
characterize effective Zeeman splitting of the interface¬ 
like subband due to in-plane and perpendicular magnetic 
fields (similar terms for the heavy hole subband are small 
enough to be neglected), and D describes spin-orbit cou¬ 
pling due to bulk inversion asymmetry. As the upper left 
and lower right 2x2 blocks of the Hamiltonian (1) cor¬ 
respond to the states with opposite projections of spin, 
the terms D and KB± are responsible for spin mixing. 

The Hamiltonian (1) can be presented in the form H = 
Hq + Hi, where 
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Since the 2D system is supposed to be a topological in¬ 
sulator, the case of inverted subband ordering, A < 0, 
is considered below. The Hamiltonian Hq describes the 
states with zero component of momentum along x axis in 
the absence of potential perturbations <p T , S r and mag¬ 
netic fields. Assume that there is a boundary at y = 0 
and the 2D system is at y > 0. The eigenstate problem 
(Hq — £)ip(y) = 0 has a couple of edge-state solutions 
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with wave functions evanescent from the boundary and 
energies inside the gap: 


entering the diagonal part of h are the edge state veloc¬ 
ity v and momentum displacement p 0 defined as 
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where it = i/(|A|/2 + e)/(| A|/2 — e) and Ak±(e) = 
\J (A/2 ) 2 — £ 2 ±i.D. To complete the eigenstate problem, 
the Hamiltonian Hq should be supplemented with bound¬ 
ary conditions which allow one to find proper orthogonal 
combinations of the wave functions (3) and correspond¬ 
ing eigenvalues of energy. Without specifying these con¬ 
ditions, let us apply a simplifying assumption that the 
particle-hole symmetry of the Hamiltonian Hq is retained 
in the presence of the boundary. This case is realized, 
for example, if the normal insulator state at y < 0 is de¬ 
scribed by a Hamiltonian of the same form as Eq. (2) but 
with large positive A; then the boundary conditions take 
the form [gj V'i(O) = — 02 ( 0 ) and ^ 3 ( 0 ) = — ^> 4 ( 0 ), where 
il>i are the components of the columnar wave function 0 . 
Under this assumption, the eigenstates are double degen¬ 
erate, with energy e = 0, so u = 1. The functions (3) 
are orthogonal to each other under these conditions. It 
is convenient to form the basis of eigenstates by using 
normalized orthogonal combinations of these functions, 
M /1 {y) and ^(y), according to 
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where n = 1,2 numbers the edge states, 

1/4 


and 
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K± = K(l±i(x), K = /X = (6) 


The quantity k is the inverse length of decay of edge-state 
wave function from the boundary. The dimensionless pa¬ 
rameter p characterizes strength of spin-orbit coupling. 

Applying this basis, i.e. searching for the solution of 
the eigenstate problem with Hamiltonian (1) in the form 
cj>i(x)' i & i(y) + <j) 2 (x)' i S? 2 (y), one obtains a 2 x 2 matrix 
equation describing ID motion in the edge states: 


{h - £)<(>(x) = 0, (7) 
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The complex energy entering the non-diagonal part is 
linear in magnetic field: 


KB X pevB z 

2 yJ \ + /i 2 1 2 + 2kc(1 + /i 2 )' 


(9) 


If the magnetic field is directed in the {xz) plane, g B 
is real. The contribution of the magnetic field into the 
effective Hamiltonian h has the general form given in Ref. 
12. However, in contrast to Ref. 12, all the parameters 
are now explicitly related to the quantities entering the 
Hamiltonian of 2D electrons in HgTe quantum wells, Eq. 

(I)- 

The usage of truncated basis H f i(?/) and M/ 2 ( 2 /) (which 
is actually the key approximation in derivation of any 
effective Hamiltonian) provides limitations for applica¬ 
tion of Eq. (7). This equation cannot be used in the 
region of energies beyond the 2D gap, |e — ip x \ > |A|/2, 
where the extended 2D states exist together with the 
edge states and can mix with them. Such mixing occurs 
even at \e — ip x \ < |A|/2 if the fluctuations of the gap, 
d r , are strong. For this reason, |<5 r | has to be small com¬ 
pared to |A|. In fact, the gap fluctuations do not enter 
Eq. (7), because the latter is obtained in the first order 
of perturbation theory in the Hamiltonian H\. Calcula¬ 
tion of the second-order contributions leads to additional 
terms associated with S r , such terms can be identified 
as coordinate-dependent corrections to the quantities k 
and p entering v, po , and g B in Eq. (7). Since these 
corrections are assumed to be small, Eq. (7) is justified. 
In spite of the limitations described, the expression (8) 
for the edge state velocity v is exact (valid for arbitrary 
energy) at B = 0. 

Equation (7) is similar to one-dimensional Dirac equa¬ 
tion for massive spinless particles. Consider the main 
properties of this equation. First of all, the terms ±i>po 
entering the diagonal part of h describe a shift of electron 
spectrum in momentum space, due to the perpendicular 
magnetic field. These terms have no effect on edge state 
propagation [9!] and can be removed from Eq. (7) by 
a substitution <j){x) —> <p(x) exp(ipoX/h) equivalent to a 
gauge transformation. For this reason, the terms ±i>po 
in h are omitted below. 

The free motion in edge states, when ip x = 0, is char¬ 
acterized by constant momentum p and described by the 
energy spectrum 


£ P = ± s/v 2 p 2 + |.9b| 2 , (10) 


where 0 = (</>i,</> 2 ) T is the columnar wave function of 
the edge states and ip x = 2 k dyip r e~ 2K,v is the elec¬ 
trostatic potential acting at the edge. The parameters 


which has a gap equal to 2( 93 1- Both the in-plane field 
and the perpendicular field contribute to the gap. How¬ 
ever, the contribution of the perpendicular field appears 
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only due to orbital effect of this field and in the pres¬ 
ence of spin-orbit coupling Q. when p, / 0, while the 
contribution of the in-plane field does not require this 
coupling M- In the case of tilted field with B y = 0 and 
B z /B x = —ncKyJ 1 + Y/evp the gap disappears. Mak¬ 
ing estimates for HgTe quantum well of width 7.3 nm 
(|A| ~ 20 meV) and assuming D ~ 1 meV, one finds 
that the ratio B z /B x corresponding to zero gap is of the 
order of unity. The origin of gap disappearance can be 
explained as a mutual compensation of the orbital and 
Zeeman effects of the magnetic held, which are caused by 
B z and B x , respectively. 

If the gap is equal to zero, Eq. (7) has two exact 
solutions: 

0 (1) ( x ) = C 0 ) ’ 0 (2) ( x ) = ( ) > ( n ) 



Since the probability current for the states of Eq. (7) 
is equal to u(|</> 2 (a;)| 2 — |</>i(x)| 2 ), these solutions de¬ 
scribe constant probability currents equal to — v (left- 
moving) and v (right-moving), respectively. Regardless 
to the form of the potential ip x , there is no backscat- 
tering. Thus, the approach presented above relates the 
topological protection of edge states to a basic property 
of massless ID Dirac model. This property is responsible 
for ideal transmission of electrons through potential bar¬ 
riers (Klein paradox jl7|), observed in graphene at nor¬ 
mal incidence [TH, Qj|. In the case of nonzero gap, the 
transmission depends on details of the potential, as illus¬ 
trated below in some examples of quantum-mechanical 
problems. 


III. PROPAGATION OF GAPPED EDGE 
STATES 


Consider propagation of electrons in the edge states 
described by Eq. (7) through finite-size regions with 
inhomogeneous electrostatic potential. Since the meth¬ 
ods of solutions of similar quantum-mechanical problems 
for Dirac model are well known [a, and the subject is 
under increasing attention in connection with studies of 
graphene [HI, the details of calculations leading to the 
results presented below are omitted. The electric current 
carried by a single edge through an inhomogeneous re¬ 
gion is given by j = (e/h) f deT E [f L (e) - / fl (e)], where 
II (fa) are the energy distribution functions of electrons 
in the edge states on the left (right) sides of the poten¬ 
tial perturbation and T e is the transmission probability 
related with the reflection probability R s in the usual 
way, T e + R e = 1. The simplest important cases are the 
following. 

Potential step. In the case of the potential step ( ip x = 0 
at x < 0 and ip x = V at x > 0) shown in Fig. 1 (a) the 
transmission of electron with energy e is nonzero if | V — 




FIG. 1: Transmission through a potential step (a), a potential 
barrier (b), interband transmission in a homogeneous electric 
field created in a Zener diode structure (c), and transmission 
through a magnetic barrier defined by a local gap at |a:| < 6/2 
(d). Solid lines plotted at ± |(?b| denote the potential 
energies of the conduction and valence bands of edge-state 
ID electrons; in the part (c) the band edges of 2D electrons, 
at ± |A|/2, are also plotted. 


e| > |5 b|• Introducing Ai = \g&\/e and A 2 = \gu\/{V-£) 
one gets 


Re 


\m + 52I 2 n _ K 
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( 12 ) 


This expression is formally valid also at \V — e\ < |<7b |, 
when 772 is a complex quantity; in this case the reflec¬ 
tion is perfect, R e = 1. The expression (12) is con¬ 
siderably simplified for a special condition e = V/2 
(symmetric n — p junction, as shown in Fig. 1 (a)): 
R e = 4|< 7 b| 2 /A^ 2 . If the gap is small so that |A»| <C 1, 
the reflection probability is small and quadratic in B: 
Re^\gB\ 2 V 2 /[2e(V-e)} 2 . 

Potential barrier. The case <p x = 0 at |x| > 6/2 and 
<fix = V at \x\ < 6/2 represents the potential barrier 
shown in Fig. 1 (b). The reflection coefficient is 


„ sin 2 (kb) (1 ~ Aj)(l — A|) 

£ sin 2 (fc6) + Q ’ (Ai + A 2 ) 2 


(13) 


where k = y/iV — e) 2 — \g&\ 2 /hv is the wavenumber of 
electron inside the barrier. If \V — e| > |5 b| (e = £1 
in Fig. 1 (b)), electrons can propagate in the barrier 
region, k is real, and R e demonstrates interference os¬ 
cillations and goes to zero each time when kb /it is inte¬ 
ger. Since k is controlled by the magnetic field, the os¬ 
cillations can be induced by sweeping B. If \V — e| < 
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|5b| (e = £2 in Fig- 1 (b)), electrons cannot prop¬ 
agate in the barrier region. In this case k is imagi¬ 
nary, k = i\k\, sin 2 (kb) = — sinh 2 (|fc|6), and transmission 
through sufficiently wide barriers is exponentially small, 
T e oc exp(— 2\k\b). If |A*| -C 1, the reflection is small, 
R e ~ \g B \ 2 V 2 sin 2 (kb)/[£(V - e)] 2 . 

Homogeneous electric field. The motion of a Dirac elec¬ 
tron in the potential tp x = Fx assumes transition of this 
electron between conduction and valence bands in the 
region near the point e = ip x . For massless Dirac Hamil¬ 
tonian this transition takes place with unity probability. 
For massive Dirac Hamiltonian the electron undergoes 
interband (Zener) tunneling, which is formally described 
within an exactly solvable problem of tunneling through 
a parabolic barrier |[2Q]. The energy-independent trans¬ 
mission probability is given by the expression 

T = exp (— n\gB\ 2 /hFv) . (14) 

A similar expression describes Zener tunneling in 
graphene at nonzero angle of incidence [21] . If the mag¬ 
netic field is weak so that |5 b| 2 HFv/tt , one gets a 
small reflection probability quadratic in magnetic field, 
R = 7T|,9b | 2 /hFv. The homogeneous electric field can be 
created, for example, in planar n—p junctions IaI formed 
in HgTe quantum well structures with split gates, Fig. 1 
(c). In the case of normal insulator the conductance of 
such a junction is proportional to L y exp(—7rA 2 /4F'A), 
where L y is the sample width. In the case of topological 
insulator, the edge state transport causes an additional 
contribution to conductance 122] . which is equal to 2 e 2 /h 
in the absence of magnetic fields and can dominate over 
the normal (bulk) contribution as the latter is exponen¬ 
tially small at A 2 FA. If the magnetic field is present, 
the edge-state contribution is multiplied by a factor of T 
from Eq. (14). Since F ~ |A|/d, where d is the length 
of the inhomogeneous region, the factor in the exponent 
of Eq. (14) is estimated as — 27 tkc%b/A| 2 . For d of 
micrometer scale, one has nd 1, because k - 1 ~ 30 
nm. This means that, despite of the assumed smallness 
of the ratio |^b 1/1A|, an exponentially small edge-state 
transmission in planar n — p junctions is feasible in the 
presence of a magnetic field. Therefore, application of 
a magnetic field can be used for efficient control of the 
current through Zener diode in the topological insulator 
state. 

Notice that Eqs. (12)-(14) coincide with correspond¬ 
ing results for one-dimensional Dirac model E2 , where 
|<7b| plays the role of a ’’mass parameter”. The reason 
for this is the following: although the effective Hamilto¬ 
nian h is formally different from one-dimensional Dirac 
Hamiltonian, it can be reduced to the latter by a unitary 
transformation. 

Magnetic barrier. The case of inhomogeneous mag¬ 
netic field is described by Eq. (7) with x-dependent <?b- 
Consider a rectangular magnetic barrier, implying that 
the field is absent at |x| > 6/2 and finite at \x\ < 6/2. 
This leads to the potential configuration where the gap 
exists locally at |x| < 6/2, as shown in Fig. 1 (d). The 


reflection probability is 


Re 


h 2 v 2 k 2 \ 
I ,9b | 2 sin 2 (kb)) 


(15) 


where k = \Je 2 — |.9b| 2 /7w. Similar to the case of po¬ 
tential barrier, electrons either can or cannot propagate 
in the barrier region: |e| > |<7 b|, real k or |e| < |<7 b|, 
imaginary k , respectively. In the first situation the trans¬ 
mission probability shows interference oscillations, while 
in the second one the transmission can be exponentially 
small. The model of magnetic barriers is suitable for de¬ 
scription of electron scattering by magnetic impurities in 
the absence of magnetic field. Since the spin-dependent 
interaction of electron with such an impurity occurs at a 
short range, it is natural to assume small 6 in Eq. (15) 
for such applications. If the product kb is also small, 
Eq. (15) describes a small energy-independent reflection 
probability R ~ | <7b | 2 6 2 //i. 2 v 2 . 

The validity of the results of this section is determined 
by the range of applicability of the effective Hamilto¬ 
nian (7). In particular, an estimate of the second-order 
perturbation contributions to h allows one to conclude 
that the results are applicable in the range of magnetic 
fields where the energies \KB\ and \K Z B Z \ are small in 
comparison to |A|. For HgTe wells both K and K z are 
estimated as 1 meV/T, so the results are justified within 
1 Tesla range provided that |A| > 10 meV. 


IV. DISCUSSION 

The opening of the gap in the edge-state spectrum of 
2D topological insulators under the action of magnetic 
field is the fundamental consequence of violation of the 
time reversal symmetry. A one-dimensional quantum- 
mechanical equation derived in this paper describes mo¬ 
tion of electrons in the edge states in these conditions. 
The contribution of the in-plane field to the gap en¬ 
ergy 2|5 b| arises because of Zeeman interaction, while 
the contribution of the perpendicular field exists purely 
due to orbital effect of magnetic field in the presence of 
spin-orbit coupling. The gap is sensitive not only to the 
strength of the magnetic field but also to the direction of 
this field, and vanishes at a special field orientation. 

When the gap is present, the scattering of electrons 
between counter-propagating edge channels can be re¬ 
sponsible for suppression of the ballistic conductance in 
HgTe quantum well samples. Let us consider a model 
of transport in a disordered HgTe quantum well near the 
charge neutrality point, when Fermi energy is in the mid¬ 
dle of the gap between 2D subbands. Since the magnetic 
field creates a small gap 2\gs\, an edge-state electron 
moving in a spatially smooth random potential tp x , see 
Fig. 2 (a), experiences multiple Zener tunneling events 
for which the averaged potential energy tilt F is esti¬ 
mated as 2w/l ci where w is the mean amplitude of the 
potential variations (w is assumed to be much larger than 
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(b) 
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FIG. 2: The energies of the band extrema for edge-state ID 
electrons, given by the expression ip x ± |(/b|, in the cases of 
smooth (a), and sharp (b) random potentials. Electrons with 
energy e experience multiple Zener tunneling (a) or multiple 
transmission across potential steps (b). In both cases the 
reflection probability is quadratic in magnetic field. 


|(?b|) and l c is the mean distance between adjacent local 
minima and maxima of the random potential. Each tun¬ 
neling event occurs with small reflection probability, but 
the total number of such transitions is large and esti¬ 
mated as N = L/l c , where L is the sample length. Ne¬ 
glecting interference effects, one may roughly estimate 
the total reflection probability as NR, with R = 1 — T 
taken from Eq. (14). Then, the relative increase in re¬ 
sistance is expected to be of the order of Ti\g-B\ 2 L/2hvw. 
If the magnetic field is perpendicular to the quantum 
well plane, this dimensionless quantity is estimated as 
2itA 3 D 2 L/\N\ a £ a w, where £ = ^hc/e\B z \ is the mag¬ 
netic length. Using A = 360 meV nm, |A| = 20 meV, 
D = 1 meV, and assuming w = 3 meV and L = 2 gm, 
the relative increase in resistance is given by 2.5[I?2(T)] 2 . 
If the random potential is sharp, as shown in Fig. 2 (b), 
the electron is transmitted across a number of potential 
steps, and the relative increase in resistance is estimated 
as \g-B\ 2 L/w 2 l c , where l c is now the average distance be¬ 
tween potential steps. Assuming l c — 10 nm, one gets 
a stronger effect of the magnetic field on the resistance 
compared to the case of smooth random potential. How¬ 
ever, the model of smooth potential seems to be a more 
realistic one. 

Therefore, according to the theory described in this 
paper, the increase in resistance is proportional to the 
square of the magnetic field and becomes essential in the 
fields of the order 1 T. Since the contributions of the 
in-plane field and the perpendicular field to the gap en¬ 


ergy are comparable (see Sec. 2), the effect of both fields 
on backscattering should be similar. In agreement with 
this theory, the numerical simulations of transport on a 
2D disordered site model [§] show a weak, quadratic in 
magnetic field increase in resistance in the case of weak 
disorder. For strong disorder, when the potential ampli¬ 
tude exceeds the bulk gap energy, the simulation gives a 
completely different result: a strong, linear in B z increase 
in resistance and much weaker effect of the in-plane field. 
Since it is the behavior actually observed experimentally, 
one may conclude that the strong disorder is a necessary 
precondition for the dramatic influence of the perpendic¬ 
ular magnetic field on ballistic resistance and that the 
strong disorder is present in the samples investigated. 
The case of strong disorder, when transitions of electrons 
between the edge and bulk states are possible, is beyond 
the range of validity of the effective Hamiltonian theory. 

On the other hand, the theory can explain the observed 
[3. fjj increase in resistance under the in-plane magnetic 
field of the order 1 T. This increase apparently begins 
with a quadratic field dependence and is sensitive to the 
gate voltage controlling position of the Fermi energy in 
the bulk gap @. Moreover, the dependence of the resis¬ 
tance on orientation of the in-plane field shows (Ref. 3, 
p. 68) that the field perpendicular to the current leads 
to a stronger effect. This is in agreement with the theory 
demonstrating that the field B y creates a larger gap than 
the field B x , see Eq. (9). The difference is controlled by 
the parameter g = 2D/|A|, which is not small in the case 
when the bulk gap |A| is of the order of 1 meV. Indeed, 
a considerable dependence of the resistance on orienta¬ 
tion of the field has been observed [3| for a sample with 
a small bulk gap. Numerical simulations ||J, carried out 
for g <C 1, show a negligible orientation dependence, in 
agreement with the theory. 

In summary, the parameters of the effective Hamilto¬ 
nian for edge states in two-dimensional topological in¬ 
sulators based on HgTe quantum wells are calculated 
analytically. The effect of magnetic fields on quantum- 
mechanical properties of these states in the presence of 
potential perturbations is studied. A feasibility of broad 
control of ballistic transport of edge-state electrons by 
magnetic fields is demonstrated. The results of the ef¬ 
fective Hamiltonian theory are compared to both exper¬ 
imental findings and numerical simulation data on mag¬ 
netoresistance of HgTe quantum wells. The theory ex¬ 
plains the observed behavior of resistance under the in¬ 
plane magnetic field, but does not explain a strong lin¬ 
ear increase in resistance as a response to perpendicular 
magnetic field. This increase is likely caused by the tran¬ 
sitions of electrons between edge and bulk states, which 
are not accounted by the effective Hamiltonian theory 
and are possible in the presence of high-amplitude disor¬ 
der potential. 
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